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We present an optomechanical accelerometer with high dynamic range, high bandwidth and read¬ 
out noise levels below 8 p.g/\/Hz . The straightforward assembly and low cost of our device make it 
a prime candidate for on-site reference calibrations and autonomous navigation. We present experi¬ 
mental data taken with a vacuum sealed, portable prototype and deduce the achieved bias stability 
and scale factor accuracy. Additionally, we present a comprehensive model of the device physics 
that we use to analyze the fundamental noise sources and accuracy limitations of such devices. 


I. INTRODUCTION 

Accelerometers and gyroscopes form the fundamental 
building blocks of inertial sensing mm- Devices with 
a wide range of bandwidth, precision, accuracy and dy¬ 
namic range are available and deployed in various appli¬ 
cations, including commercial products, medical devices 
[3], construction engineering [3], natural resource explo¬ 
ration [5], inertial sensing for autonomous navigation [6] 
and fundamental research [71IH]. The critical parameters 
of an accelerometer with a given bandwidth are its pre¬ 
cision, often denoted as acceleration noise Oacc; its bias 
stability, the long-term drift of the DG acceleration absl 
and the scale factor S uncertainty, which describes errors 
in the value relating physical observable quantities, such 
as the output voltage Uacc to acceleration {va^cc = S • Oext, 
(S) = 1 V/m/s^ ). 

Gommercially available accelerometers require labo¬ 
rious calibrations to ensure a minimal deviation of 
the above described parameters. Such calibrations 
are often performed at National Metrology Institutes 
(NMIs) [9], where the test unit is mounted onto an 
interferometrically-interrogated reference shaker system 
[inHE]- Such calibrations reach relative uncertainties 
of the order of 10“^ — 10“^'^. Already calibrated de¬ 
vices, known as reference accelerometers, can, in turn, 
be used to calibrate further devices using more simple 
back-to-back measurement set-ups. However, such de¬ 
pendent calibrations are accompanied by an inevitable 
degradation in uncertainty. 

Optomechanical accelerometers interrogated using 
fiber interferometric methods have recently demon¬ 
strated high levels of readout precision [imn], using 
optical instead of electro-static readout. The devices in¬ 
vestigated by Guzman et al. m use monolithic fused sil¬ 
ica in-plane oscillators interrogated by fiber optic micro¬ 
cavities mni]. These devices have achieved levels of ac¬ 
celeration measurement sensitivities below I00ng/-\/Hz 
and, more importantly in this context, they provide di- 
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rect traceability to SI units, a property we denote as 
“self-calibrating”. This is achieved in a two-step pro¬ 
cess: First, external accelerations of the oscillator are 
converted to displacement via its transfer function, char¬ 
acterized by its resonance frequency ojq and mechanical 
quality factor Q. Second, displacement is turned into a 
readout voltage via laser interferometry. Similar devices 
have also recently been used as optomechanical force sen¬ 
sors for atomic force microscopy (AFM) |I8) . 

In this article, we present a similar device to Ref. [15] 
that is interrogated by a low-finesse cavity formed by a 
gap of approx. 50 |xm between two flat-cleaved fibers [Sl¬ 
its simple and cost-effective construction, combined with 
self-calibration, high dynamic range, and high bandwidth 
(exceeding 10 kHz) make this device a promising candi¬ 
date to perform autonomous navigation and on-site accel¬ 
eration calibrations of other accelerometers. We demon¬ 
strate the portability and accuracy of such a device by 
vacuum packaging it and by characterizing its readout 
noise and bias stability. Without additional feedback, 
the low finesse readout allows our device to have a scale 
factor uncertainty lower than 10“^ for accelerations from 
mg to those exceeding 2g. We also present our current 
understanding of the device physics, optics and mechan¬ 
ics and discuss the expected limits for acceleration read¬ 
out noise, accuracy and self-calibration and, in addition, 
we extrapolate this analysis to a high-finesse cavity read¬ 
out, paving the path for substantial improvements. 

II. ACCELEROMETER PROTOTYPE 

Our accelerometer device (shown in Figure uses a 
similar in-plane, monolithic fused-silica oscillator design 
described by Guzman et. al. m- It is frit-bonded onto an 
additional micro-machined fused-silica part that contains 
a relief for enabling free test mass motion. This part is, 
in turn, frit-bonded onto a 4 cm radius quartz plate with 
a thickness of 3 mm. The accelerometer is enclosed by 
a glass bell, which is glued with Torr Seal onto the base 
plate. Two flat-cleaved fibers are fed into the bell at a 
cut-out and their ends are glued into v-grooves next to 
the test mass. Flat-cleaved fibers glued onto opposite fac¬ 
ing v-grooves on the test mass act as cavity end mirrors. 
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FIG. 1: (a): Sketch of the acceleration sensing using the me¬ 
chanical oscillator and a single low-hnesse fiber micro-cavity, 
(b): Photograph of the vacuum sealed accelerometer proto¬ 
type with two hber interferometers, one mounted in the lower 
v-grooves for redundancy. 


The acceleration-driven test mass motion translates to 
optical phase changes of the fiber micro-cavity, which is 
read out in reflection by monitoring the reflected power. 
The fiber ends acting as cavity mirrors, have a low reflec¬ 
tivity of approximately 4 %. The oscillator was designed 
to have a mass of 25 mg and a resonance frequency of 
about 10 kHz. 

After assembling the fiber cavities, we vacuum sealed 
our device by pumping on a glass pipe originally con¬ 
nected to the top of the bell and performing a vacuum 
pulling by flame heating the pipe. This is evident by the 
residual glass structure on top of the bell shown in Figure 
[^. The device is now available as a portable accelerome¬ 
ter which we connected to our measurement set-up shown 
in Figure 



FIG. 2: Experimental set-up for the accelerometer charac¬ 
terization. A laser beam from a widely tunable laser m is 
send into one of the fiber cavities and the reflected power is 
detected on a photo receiver {vr). To increase optical power 
for the performance test an erbium-doped fiber amplifier |20) 
(not shown) was integrated after the laser source. Part of 
the laser light is split off to stabilize the laser power by ac¬ 
tuating on a fiber-based electro-optic amplitude modulator 
(EOAM). A fiber-based unequal arm-length Mach-Zehnder 
interferometer, with 10 m hber delay between the arms, is 
used as frequency reference. It is read out via a balanced 
detection scheme and the error signal is used to stabilize the 
laser frequency by actuating on the laser pump current and on 
the laser cavity piezo. The accelerometer is placed on an iso¬ 
lation platform to reduce coupling of unwanted accelerations. 
A piezo crystal (PZT) mounted on the platform is used to ex¬ 
cite the accelerometer. The second, unused output port of the 
50/50 beam splitter was additionally terminated by wrapping 
the hber with a very small bending radius and thereby intro¬ 
ducing excess losses of the light traveling through the hber 
core. 


A. Laser wavelength scan 

We use a widely tunable laser m to measure the de¬ 
pendency of the reflected light power Pr, and the corre¬ 
sponding photo detector output voltage vr at the laser 
wavelength A. Figure|^shows the response for one of the 
two cavities, which was used throughout this study. 

We first fit our data using the following sinusoidal 
model approximation of the fiber cavity response. 

Vr^s = VdG - VDC ■ kcos{4:Tt^), ( 1 ) 

where z;dc is the DC output voltage measured in reflec¬ 
tion, K is the optical contrast, and Zm is the macroscopic 
cavity length. By fitting our measurement data with 
this response we determine vdc = 2.85 V, k = 0.97 and 
Zm ~ 42.6 pm. We found a low uncertainty of about 
0.15% for the slope at the critical quadrature points of 
the response between the fit and the data. The residuals 
indicate that our response function models are not able to 
fully describe the measured behavior. We attribute this 
to effects related to diffraction, residual misalignments 
and coupling efficiency into the fiber mode |21j . which 
become relevant due to our short, non mode-matched 
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FIG. 3: The top shows the measured cavity response vr to a 
wavelength sweep in comparison to the fitted sinusoidal (wi?,s) 
and Airy function (vR,a) models. The bottom plot shows the 
residuals of the two fits. The gray, vertical lines indicate the 
range of data used for the fits. The laser power was actively 
stabilized during the sweep to account for a changing laser 
output power. The inset in the top right corner shows a zoom 
into the quadrature point for the response and the residuals. 
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FIG. 4: Oscillator ring downs in comparison to the fitted 
envelope functions. The measured data was high-pass filtered 
to exclude low frequency noise and to make the decay clearly 
visible. 


cavity and the wide tuning of laser wavelength. However, 
the impact of these effects is much lower around our op¬ 
erating point. This is discussed some more in Section [V| 
By tuning the laser wavelength to a quadrature point of 
the response (Ag), we maximize the ratio between voltage 
and length change to a now calibrated value [55] . 

47r 

l^VR = t'DC • K — AZm- (2) 

Ag 

We also analyze our measurement data with the Airy- 
function, which fully takes multiple reflections into ac¬ 
count [TT] : 


1 - cos {^Zm 
l-hi?2_2i?cos(^ 

^(3) 

Here we define an offset voltage Uoff, an amplitude scal¬ 
ing voltage v~f that includes the optical contrast and the 
reflectivity R « 0.04. Around the now slightly shifted 
point of maximum slope, Aq, we can again define an ef¬ 
fective linear coupling of displacement into voltage in a 
more general form. 



^R,a — ^off 4“ ^7 


{i + nr 


Azm- (4) 

Zm \dXJ 

B. Ring down 

To estimate the oscillator mechanical characteristic pa¬ 
rameters we perform a ring-down experiment with our 
device. We excite it close to the resonance frequency, 
turn off the excitation, and then monitor the oscillation 
decay. A measurement for our device is shown in Figure 
We fit our data using the formula 

_ UJQt 

vrd = Vo ■ e~^ sin {ujot + (j)). (5) 


From this fit we estimate the natural frequency of the 
oscillator ojq (« 27r x 10646.56 Hz) with a fit standard 
error on the order of 10“®. We then perform a software 
demodulation of the signal at ujq to extract only the ex¬ 
ponential decay of the oscillation amplitude Vrd, which is 
governed by the following equation. 

_ 

Vrd = Vo ■ e . (6) 


By fitting this formula to the decay we determine the 
quality factor Q « 12379.07 with a standard uncertainty 
of about 15 X 10“®. The resulting fit curve is shown in 
Figure]^ as the envelope of the ring down. A correction 
of Wo due to a damping-induced frequency shift is not 
necessary, since the correction factor is fully negligible at 
these levels of Q. To further refine and test the long term 
stability of wq and Q we perform the above described 
analysis for multiple, consecutive ring down measure¬ 
ments. For 250 ring downs, spaced over 25 minutes (1500 
seconds) of measurement time, we determine a standard 
deviation for wq on the order of 0.1 x 10“® and 200 x 10“® 
for Q, each an order of magnitude larger than the single¬ 
shot statistical variation. Using the transfer function of 
a damped harmonic oscillator, we have the relation be¬ 
tween external accelerations and displacement: 


THo(<j-’) 


Az(w) 

^ext (^) 


1 

Wg — W^ -b i^W 


( 7 ) 


C. Noise performance 

With our device resting on a vibration isolation plat¬ 
form we measure the spectra of the output voltage to 
determine the readout noise floor. Using equations]^ and 
[^we convert the voltage spectra into the corresponding 
displacement. The results of this are shown in Figure 
[^for different laser stabilization schemes (see Figure]^. 
With both laser amplitude and frequency stabilization 
the spectrum shows the thermally excited peak of the 
oscillator resonance, which sticks out of an almost flat 
noise floor of about 15fm/-\/Hz . The noise increases at 
lower frequencies, though the magnitude is strongly re¬ 
duced by the stabilizations. We estimate the shot noise as 
explained in Section [HI Fj with a measured DC power on 
the photodiode of 60 |xW. The dark noise measurement 
includes contributions from the photo receiver and the 
measurement devices. This was measured by blocking 
all light going onto the measurement photodiode. The 
laser frequency stabilization allowed us to strongly reduce 
the excess frequency noise, likely induced by coupling 
of acoustics and vibrations into the tunable laser head 
[^|5Tj. The unequal armlength Mach-Zehnder config¬ 
uration is able to operate at any given wavelength and 
can, therefore, be used with the tunable laser, necessary 
to perform the self-calibration. 

Using the harmonic oscillator transfer function model 
of the mechanical oscillator we can determine the corre- 
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FIG. 5: Measured displacement spectra without stabiliza¬ 
tions (R), with amplitude stabilization (i?a), frequency stabi¬ 
lization (i?f) and both stabilizations (i?f&a) running. Shown 
are also the estimated shot noise level {Rn) and the measured 
dark noise (i?dn). The thermally excited resonance peak of 
the mechanical oscillator is visible at /q. 



FIG. 6: Measured acceleration spectra with (/?f&a) and with¬ 
out stabilizations (R) running. Shown are also the estimated 
shot noise level (Rn) and the measured dark noise (i?dn). 

spending acceleration noise spectra for a given displace¬ 
ment noise. Using a data acquisition system we per¬ 
formed long-term measurements. We converted the mea¬ 
sured voltage into acceleration by filtering and scaling 
it accordingly. The resulting acceleration noise spectra 
are shown in Figure Our device achieves an accel¬ 
eration noise floor better than 8|J.g/-\/Hz above 1 kHz. 
Its low frequency performance is limited by a 1// noise, 
leading to levels of 900 |J.g/-\/Hz at 1 Hz and better than 
40mg/-\/Hz at 10 mHz (corresponding to a displacement 
noise of 2 pm/'s/Hz and better than 90 pm/'s/Hz respec¬ 
tively) . 

D. Allan variance 

To investigate the bias stability we determined the Al¬ 
lan deviation of our readout. The results are shown in 


Figure with and without stabilization. At short inte¬ 
gration times we are limited only by shot noise. The 
exact nature of the increase at longer integration times 
is currently under investigation, and it is presumed to 
be caused by residual frequency noise or parasitic stray 
beams in our fiber set-up. Each of these contributions 
is susceptible to acoustic and thermal fluctuations, caus¬ 
ing the noise in our set-up to not be constant. This is 
evident from the second, shorter measurement with sta¬ 
bilizations shown in Figure which presents the best 
levels of long-term stability achieved during our mea¬ 
surement campaign, even though the readout noise floor 
was signihcantly higher. Our frequency reference Mach- 
Zehnder interferometer was placed in the same thermal 
environment as our laser source and was isolated pas¬ 
sively against thermal fluctuations. The use of active 
thermal stabilizations for the fiber interferometer may 
help to improve the performance in future implementa¬ 
tions. We reach the minimum of our Allan deviation at 
about 1 ms of integration time with a value of 3 x 10“^ g. 

During our investigation we found that the use of fiber 
circulators, to feed light to the accelerometer and to de¬ 
tect the reflection, caused an increased coupling of laser 
frequency noise into the measurement. We attribute this 
to parasitic beams, introduced by the excess leakage in 
the return path of the circulator, that are phase modu¬ 
lated relative to our signal of interest by the frequency 
noise. These beams, or stray light, can reach significant 
amplitudes in comparison to the reflected signal of the 
low hnesse cavity, which is only about 4% of the inci¬ 
dent power, and contaminate the interfered signal, lead¬ 
ing to enhanced coupling of laser frequency noise and 
other phase noise in our fiber set-up. Using a 50/50 
fiber beam splitter we were able to actively tune and ulti¬ 
mately minimize this coupling by terminating the open, 
unused output of the splitter. This improvement in sta¬ 
bility comes at the cost of reducing the detected power to 
about one quarter of the levels achieved with a circula¬ 
tor. Future implementations have to take these parastic 
signals into account to reduce the coupling of laser fre¬ 
quency noise to the signal from the actual cavity. 

As an example of what may be possible with such a 
device having a higher-finesse optical readout, Figure 
shows Allan deviations for the high hnesse cavity read¬ 
out implemented by Guzman et al. [15j . which were 
measured with a hnesse of 1600 and a hxed-frequency 
laser. The extended curve shows levels corresponding to 
an extension of the measured spectral density down to 
1 mHz (data was available down to 20 Hz) with a read¬ 
out dominated by an assumed laser frequency noise with 
191kHz/-\/Hz at IHz, levels that are regularly achieved 
with rather simple laser stabilization schemes, and a 
cavity length of 172 pm. The strongly reduced readout 
noise hoor for this system leads to Allan deviations below 
10-5 g. 



















5 


10" 


10“ 


lO-’ 


^ 10“" 


10 


10" 



10"^ 10“= 10“^ 10“' io“ io' 10^ 10^ 


r /sec 


The high stiffness and the correspondingly small test 
mass displacements of our device lets it be described as a 
damped harmonic oscillator (HO) with an angular reso¬ 
nance frequency cuo, a mechanical quality factor Q and an 
effective oscillator mass m. The conversion of accelera¬ 
tion into displacement is given, in the frequency domain, 
by the transfer function Tho described in Equation 
Correspondingly, a measured displacement in the time 
domain is converted back into an acceleration by com¬ 
puting 


aacc(t) = ( 8 ) 


FIG. 7: Measured Allan deviations of our detector out¬ 

put signal converted into acceleration. Shown are the results 
with both stabilizations (Ttf&a), the estimated shot noise level 
[Rn) and the measured dark noise (i?dn). A second measure¬ 
ment with both stabilizations and a higher readout dark noise 
floor (i?f&a, 2 ) is shown as well. It shows an improved stabil¬ 
ity at longer integration times, indicating that the noise is 
not constant at longer integration times. For comparison we 
also show the Allan deviations calculated from the measured 
power spectral densities for the high-finesse cavity readout 
(i?p=1600 ) m- We also extended the corresponding spectrum 
with a 1/f noise model expected for laser frequency noise of 
about 200kHz/\/Hz at IHz {R F=i600,ext) and with a cavity 
length of 100 pm (see Equation 141. 


B. Thermal noise 

The acceleration readout is fundamentally limited by 
thermal acceleration noise Uth induced by the finite tem¬ 
perature of the oscillator and mechanical losses [25]. 




-ikBTuJo 

mQ 


( 9 ) 


This is approximately a white noise that is simply added 
to Oext, and therefore, only contributes to the readout 
noise floor, but it does not influence the scale factor, nor 
the bias stability. 



FIG. 8: Simplified model of the metrology chain from exter¬ 
nal oscillator acceleration to test mass displacement. 


III. METROLOGICAL ANALYSIS 

Any external acceleration measured with our devices is 
converted into a readout voltage. In the following, we de¬ 
scribe and model this conversion and its self-calibration 
as a chain of individual transducer steps. This allows us 
to i) determine readout noise of the devices, scale fac¬ 
tor uncertainty, bias stability and dynamic range for a 
given set of parameters, and ii) to determine the param¬ 
eters that are necessary to reach a desired measurement 
uncertainty. 

A. Mechanical oscillator 

Our devices use an in-plane monolithic oscillator made 
of fused-silica to convert external accelerations Oext into 
measured test mass displacement Az (recall Figure [^. 
A block diagram of this conversion is shown in Figure 


C. Conversion into displacement 

The scaling of acceleration into displacement depends 
only on two parameters, ujq and Q, and on the signal 
frequency w, as indicated by Equation]^ We can esti¬ 
mate the scaling at two critical points of the frequency 
response. At DC and at low frequencies the scaling sim¬ 
plifies to Tho(w ^ Wo) = — I/wq, with no Q depen¬ 
dency. Near resonance the influence of Q is maximal and 
'7ho(w = Wo) = -Q/w§. 

The estimation of both parameters is critical for the 
scale factor accuracy. The measurement of wq is done 
by comparing the oscillation to a frequency standard, 
which, in turn, makes the uncertainty of the standard one 
of the parameters fundamentally limiting the achievable 
scale factor uncertainty. The measurement uncertainty 
of Qi which can be determined with different techniques, 
is only limited by the integration time and the readout 
noise floor. Larger values of Q require in general a longer 
measurement time, due to the increase in relaxation time 

r = 

D. Deviations from the harmonic oscillator 

As briefly mentioned above, the harmonic oscillator 
model is only an approximation of the real behavior of 
the sensor. A full model that allows us to estimate the 
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FIG. 9: Shown are the ground mode and the three next 
higher order modes of the oscillator, calculated using Au¬ 
todesk Inventor finite element analysis. The colour encoding 
represents the critical displacement in z-direction, the scaling 
is arbitrary. The detection area is sketched with black stripes. 



FIG. 10: Model of the metrology chain from test mass dis¬ 
placement to reflected power. 


damping, can result in a different transfer function from 
acceleration into displacement. One direct way to in¬ 
vestigate such effects is to excite the accelerometer with 
known amplitudes at various frequencies to verify the 
transfer function, an experiment that will be conducted 
at a later point. 


Linearity 


error of our approximation is beyond the scope of this 
article. However, we present the relevant parasitic ef¬ 
fects, included as additional transfer function Tpar in our 
model, by separating them into three categories and we 
discuss their expected behavior. 

Higher order modes 

The parallelogram design of our oscillator ensures that 
higher order modes are well separated in frequency and 
that their main axis of motion is perpendicular to the 
fundamental mode. Figure shows the results of a finite 
element model analysis of our sensors, which verifies the 
mode spacing, as well as the reduced coupling of higher 
order modes into the critical z-axis. As we can see from 
the plot, the detection of the oscillator motion close to the 
top surface is not ideal for reducing coupling of the higher 
modes, especially a; 2 , into z. Hence, future devices might 
use larger v-grooves or otherwise optimized geometries 
to detect the displacement at the optimal point, if these 
modes prove to be limiting in the future. 

Non-viscous damping 

Our oscillator devices are designed to be operated in 
vacuum and under small effective displacements, signifi¬ 
cantly less than 1 nm. Hence, residual gas damping and 
mechanical losses inside the oscillator are the dominant 
contributions, which are well described by a linear viscous 
damping model. Variations in these parameters are not 
at all critical, since they are incorporated into the mea¬ 
sured Q (providing that these parameters do not change 
significantly between the measurement time and the self¬ 
calibration). Non-viscous damping effects, like structural 


Simple beam deflection theory, which is part of the ba¬ 
sis for the linear harmonic oscillator model, breaks down 
for large displacements, at which point additional higher 
order dependencies on the input frequency become rele¬ 
vant and the fundamental mode loses energy, effectively 
decreasing Q. These effects limit the dynamic range of 
the oscillator, which relates, for a given readout uncer¬ 
tainty, to a maximum acceleration. One can characterize 
these behaviors experimentally, by operating the device 
under the desired maximum acceleration. For low fre¬ 
quencies such a measurement has to determine whether 
higher harmonics of the excitation frequency are present 
in the readout. For frequencies close to the resonance this 
can be combined with a measurement of Q under varying 
levels of excitation. The total deflection for accelerations 
of up to 2 g is on the order of 5 nm. This is significantly 
smaller than the thickness of the flexure in z-direction 
(255 |xm), which leads us to expect this contribution to 
be very small. 


E. Fiber cavity readout 

Our model for the propagation of the test mass dis¬ 
placement Az into reflected optical power Pn is shown 
in Figure 

Cavity length 

The total cavity length is a combination of a constant 
term zq, the test mass displacement Az, and any un¬ 
wanted, parasitic influences Zpar caused, for example, by 
thermal expansion of the sensor. Using the laser wave¬ 
length sweep measurement, the total cavity length can be 
determined in situ. Time dependent parasitic changes 
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will cause an unwanted proportional output signal, in¬ 
ducing noise and decreasing the bias stability, and, for 
significant changes, a scale factor change. 


Cavity response 

The effective cavity reflectivity for a low-finesse exter¬ 
nal fiber micro-cavity has been studied in great detail 
[2I1E61E7]. Models are available that include the influ¬ 
ence of angular misalignment and multiple reflections. In 
case of low reflectivity one can approximate this function 
with the sinusoidal response for a simple two beam in¬ 
terference, neglecting all geometric and Gaussian beam 
effects. The resulting reflected intensity R of the cavity 
is then given as 

R = Rbc ■ ^1 — K • cos(47r^)^ . (10) 

Experimentally, this response is implicitly characterized 
by the laser wavelength scan described in Section |II A| 
Since the real response of the cavity is somewhat more 
complicated, it is expected that this approximation will 
limit the scale factor accuracy. In the following, we will 
present a few equations for the sinusoidal model, which 
are simple and useful for calculating errors due to shot 
and laser intensity noise for the low-finesse cavity, and 
we will present the more general formulas for the Airy 
function model, which is applied to determine the scale 
factor and for cavities with higher finesse. For the lin¬ 
earized sinusoidal model the reflected intensity around 
the quadrature point can be estimated 

R = Ruc + AR = i?Dc + (11) 

From this approximation it becomes clear that the re¬ 
sponse depends on the laser wavelength, as well as its 
stability. Changes in A will move the operating point 
that was used during measurement characterization, in¬ 
dependent of the applied response model, which changes 
the effective scale factor and bias. 

Using the Airy function model (See Equation]^, which 
includes multiple reflections, we generalize the cavity re¬ 
flected intensity around a given operating wavelength Aq 
as 

R = RiZm, Ao) + ^ (^) ^Zm- (12) 

= i?DC + |^(^) Az™ (13) 

In case of a hypothetical perfect cavity model, the laser 
wavelength uncertainty represents the fundamental limit 
for the displacement readout accuracy. 


Nonlinearity 


For large displacements at a given uncertainty the lin¬ 
ear approximations of Equations 10 and 13 introduce er¬ 


rors in the scale factor and generate signals at higher 
harmonics of the input signal. The dynamic range of 
our device is relatively large compared to systems using 
higher values of finesse, and rough error estimates can 
be easily calculated from the deviations of the linear ap¬ 
proximations relative to the expected response functions. 
We determine a scale factor error of less than 10“^ for 
accelerations exceeding 2g. 


Laser frequency noise 


Changes in the laser frequency (which are inverse to 
wavelength fluctuations) couple into the readout as phase 
noise, causing an effective displacement noise. This is 
caused by the interference between the reflections, which 
are delayed relative to each other by twice the cavity 
length. The effective coupling for small delays in com¬ 
parison to the readout frequencies is given as 

Zf = Zm^J- (14) 

Here c is the speed of light and / is the laser frequency 
noise. 


1. Input power 

The reflectivity of the cavity is sensed with the power 
sent into the fiber. 


Pr = PRdg + par (15) 


We determine the change in reflected power for the Airy 
function model as 


APr = 



Ao 



(16) 


Intensity fluctuations 

Relative amplitude/intensity noise (RIN) has to be 
taken into account as well. Its influence can be described 
as two separate effects. The first effect is that RIN di¬ 
rectly influences the scaling factor, as evident from Equa¬ 
tion The second effect is a coupling into the readout 
as additive amplitude noise. To quantify this effect we 
write the optical power P in terms of a constant compo¬ 
nent and an additive power noise P, 


P = Po- RIN = Po + P. (17) 

The DC component of the optical signal can now be 
rewritten as a constant term and a fluctuation noise term. 

P ■ RdC = Po ■ RdC + P ■ RdC = Po,DG + Prin- (18) 








The effective length noise zp due to RIN can be esti¬ 
mated, by computing the ratio of Prin and APr/z^. 
For the sinusoidal model this corresponds to 


Zp 


PA, 

PAttk 


= RIN- 


(19) 


For the more general cavity response we can write this 
as 


z-p = 


P R(^Zm-i -^o) 


Ar 


^ (^) 


= RIN 


Aq) 




An 


( 20 ) 


Here we omit any power losses in our fiber set-up, which 
simply scale the effective power P. Experimentally, one 
can simply measure the reflected power at quadrature 
Pdc to determine the correct scaling. 


TABLE I: List of parameter values used for calculating the 
error and accuracy limits. 


parameter 

value 

unit 

LOQ 

27r X 10646.56 

Hz 

m 

25 

mg 

rpD 

0.95 

A/W 

Ao 

1558.5 

nm 

K, 

0.97 

- 

Zm 

42 

fj-in 

T 

300 

kelvin (K) 

g 

9.81 

m/s^ 

ks 

1.3806488 X 10“^^ 

m^ kgs-^K-i 

<?e 

1.602 X 10~^® 

C 

c 

299792458 

m/s 


F. Light to voltage conversion 
Photodiode 

A photo detector is used to convert the reflected optical 
signal into a photo current, via the photodiode responsi- 
tivity rpD ((?’pd) = ampere/watt = A/W). The current 
response for the Airy function model gives 

Azr = rpppF Az„. (21) 

Zm \d\ J 


Shot noise 


The DC power on the photodiode generates a shot 
noise, which can be modeled as an effective white photo¬ 
diode current noise in- This noise depends only on the 
DC input power and the photodiode responsitivity. 

in = \/2(J'ePDC’"PD (22) 


The shot noise induced displacement noise depends on 
the ratio of this current noise to Ainjzm. For the sinu¬ 
soidal model this gives 


_ 

Ain 




rpDK^STT^PDC 


(23) 


For the general cavity response this corresponds to 


Zn = 


‘Zq^R^Zni: Aq) Zn 


rpDP 


(^)Ao 


(24) 


Trans-impedance amplifier 

The conversion into readout voltage is done using a 
trans-impedance amplifier (TIA), which is often com¬ 
bined with the photodiode in a photo receiver. At this 


point, we assume that the TIA has either sufficient band¬ 
width to minimize any frequency dependent scaling ef¬ 
fects, or that the back-end corrects for this. The TIA 
is then simply characterized by a gain 7?.tia that deter¬ 
mines the current-to-voltage ratio. 


The photo receiver elements can introduce additional 
electronic noise that can spoil the measurement perfor¬ 
mance. Care should be taken to reduce any such influ¬ 
ence to negligible levels. A simple method to do this is to 
calculate the equivalent input current noise of all contri¬ 
butions in the photo receiver |28j and to compare them 
to the expected shot noise level (see Equation 22). 


IV. ERROR BUDGET 


The overall error budget for the acceleration readout 
performance is compiled for the parameters listed in Ta¬ 
ble U 

The parameters are either chosen by design or they 
represent experimentally reproducible values. For the 
major error contributions we can now derive the nec¬ 
essary values of the other readout parameters to achieve 
a sensitivity of 10 pg/vldz . The results are summa¬ 
rized in Table [n] Thermal noise is completely negligible. 
Shot noise and RIN are broadband and can both domi¬ 
nate the high frequency behavior. The influence of RIN 
can, however, be easily reduced by implementing an am¬ 
plitude stabilization. Laser frequency noise is expected 
to dominate at low frequencies for free running lasers, 
as observed. For a well-stabilized laser source the low 
frequency performance could at some point also be dom¬ 
inated by parasitic displacement noise, induced, for ex¬ 
ample, by thermal fluctuations that drive the non-zero 
coefficient of thermal-expansion of the device. 
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TABLE II: Relevant noise sources and boundary conditions 
necessary to achieve a sensitivity of better than 10 vi.g/\/Hz . 

error value for reaching o 

source symbol parameter < 10 p.g/-\/Hz 

thermal Uth Q > 0.005 

displacement Ud Zd < 22fm/\/Hz 

frequency a/ /„ <101kHz/\/Hz 

amplitude ap RIN < 1.7 x 10 “^/'\/Hz 



TABLE III: Overview of the components that limit the fun¬ 
damental scale-factor accuracy 


transducer 


mech. oscillator 


interferometer 


parameter 


tuo, Q 


dt; \ 


SI reference 


/ref 


Ao 


simple model (viscously) damped 
harmonic oscillator 
Equation 


two-beam interference 
(sinusoidal mod el) 
Equation 10 


extensions 


higher order modes 
non-viscous damping 
non-linearity 


multiple reflections 
fiber coupling efflciency 
non-linearity 


stability 

influence 


thermal expansion 
clamping 


thermal expansion 
intensity noise 


FIG. 11: Model of the metrology chain from external to mea¬ 
sured acceleration that is used to determine the accelerometer 
performance. 


external acceleration, they determine the total readout 
error. 


^err — ftext ‘ “b Oi- 


(29) 


V. SELF CALIBRATION & ACCURACY 

We evaluate the performance of our accelerometer and 
its self calibration based on the model shown in Figure 
[m Any external acceleration experienced by our device 
is contaminated by some measurement noise a, which we 
discussed in the previous section, and is converted into 
a voltage with the actual physical scale factor, or, to be 
more precise, with the scale operator Sp. A constant bias 
Ub is added to form the total measured output voltage v^. 
During post processing we subtract an estimated voltage 
bias Ub,e and we convert into measured acceleration Om 
by applying the inverse of the estimated scale operator 
In this model we neglect any influence of scale 
factor inaccuracies on the determined noise /error levels. 
We can now write down the measured acceleration. 


— tSe ("Cb Ub,e “b ^p(c^ext “b tt)) 


(25) 


We can write down the full analytical expression for the 
scale operator, based on Equations and 




du \ 


-1 


Ao \dXj 


dt^ 


//I A 

Q dt 


(26) 


Assuming we can subtract DC biases well and neglect¬ 
ing the influence of the scale factor inaccuracies on the 
noise/error levels, we simplify Equation 25 to 

^in ~ ^p^ext “b Ci. (^/) 

If we now compute the effective error of our acceleration 
measurement (oerr = Om — Oext) we get 


— Oext • (^e 5p — 1) -b 


(28) 


We call Sg — I the scale factor accuracy 5^ operator, 
and, together with the noise and the dynamics of the 


A. Scale factor accuracy 

To discuss the estimated scale factor accuracy achiev¬ 
able with our device we determine the combined stan¬ 
dard uncertainty Ug for S^- We summarize the relevant 
effects, models and fundamental parameters in Table |lITj 
split by the two characterization steps that we use to 
estimate For a given model of the behavior of 

our system we can calculate the scale factor uncertainty 
based on measured, or estimated uncertainties of the rel¬ 
evant parameters. For a number of N parameters x with 
relative uncertainties Ux and estimated values Xg, given 


as Xg = x{l Ex)^ applied with Sg 
combined standard uncertainty as 


, we estimate the 


Ug = 




(30) 


To give an example, we evaluate the coupling of wq into 
the standard uncertainty for low frequency signals. For 
this case simplifies to 


-b 1 — 


^0,e 

, ,2 

iUn 


— (1 + £wo)^ — 1 


■ 2ea;o + ■ 


(31) 


Computing the derivative and omitting terms of the order 
we determine a scaling factor of 


dS„ 

d£u„ 


(32) 


We have summarized the uncertainties of the relevant 
parameters and their coupling factors in Table |IV| For 
the mechanical oscillator part we include two cases, at 
low frequencies and at resonance. For the interferometer 
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we include the uncertainties for the here-presented low- 
finesse readout and, for comparison, the values for the 
high-finesse readout from the earlier study m- In the 
following we discuss the individual contributions. 


Mechanical oscillator 

Fundamentally, the scale factor is limited by the vari¬ 
ations of ujq and Q, due to temperature changes, causing 
thermal expansion or material property changes, and due 
to changes in the clamping or mounting of the accelerom¬ 
eter, which also influence Q. From our earlier spec¬ 
troscopy measurements we find that ujq is stable within 
1 X 10“® in our laboratory thermal environment and Q 
can be estimated with an uncertainty of about 200x 10“®. 
In future studies we will investigate the temperature in¬ 
fluence further, to get a more detailed understanding of 
the thermal environment necessary to achieve even better 
scale factor accuracies. 


Interferometer 

The calibration of the optical metrology chain com¬ 
bines many of the earlier described transducer steps into 
a direct response of output voltage over laser wavelength, 
which is used to determine the effective ratio of volt¬ 
age over displacement The fundamental limitation 
of this approach is given by the wavelength uncertainty. 
The currently used laser system provides a resolution of 
0.01 nm during a wavelength scan. From this we can esti¬ 
mate the wave length uncertainty limit to ua = 7 x 10“®. 
Changes of the cavity length due to thermal expansion 
or of the reflected power, due to intensity noise, influ¬ 
ence the scale factor over time. Power stabilizations that 
reach levels of better than 10“® over long time scales are 
feasible and can be used, if higher scale factor precisions 
are aimed for. 

We have characterized the current experimental lim¬ 
its for our cavity readout characterization using the Airy 
function model described in Equation]^ and the resulting 
^ (see Equation 1^. Within this model the uncertain¬ 
ties of four parameters are relevant: u\, ur and 

. The standard errors of our Airy function fit point 
to an uncertainty for Zm of better than 10 x 10“®. Tak¬ 
ing the fundamental limit of u\ for this determination 
into account, we estimate uzm « 12 x 10“®. The re¬ 
flectivity of our fiber ends is only poorly estimated by 
the fit to Ur « 10%, but we can determine it better us¬ 
ing additional reflection and transmission measurements 
to about Ur « 0.5%. The coupling of reflection vari¬ 
ations into the scale factor is complex and is given by 
the derivative of the Airy function we estimate the 
uncertainty scaling to « 0.003. However, our dominat¬ 
ing noise term arises from uncertainties in estimating u.y, 
which, based on our fits, can be determined with an un¬ 
certainty of Uy^ « 1.5 x 10“^. We can exclude an influ- 


TABLE IV: List of the individual relative uncertainties and 
their scaling, together with the derived combined standard 
uncertainties for signals at low frequencies and at oscillator 
resonance. 


uncertainty 

harmonic oscillator 



w <C Wo 

iJ = UjQ 



value 

(JS rr 

value 

de-T 

Uljq 

0.1 X 10“® 2 

0.1 X 10“® 

1 

UQ 

< 0.2 X 10"^ 0 

< 0.2 X 10"^ 

1 



uncertainty 

Fabry-Perot interferometer 



R = 4% 

F = 1600 



value 

value 

dS^ 

Ux 

7 X 10“® 1 

7 X 10“® 

1 


5 X 10“® 1 

7 X 10“® 

1 

Ur 

5 X 10“^ 0.003 

78.4 X 10"® 

1.33 


1.5 X 10“® 1 

309 X 10"® 

1 



Uc 

R = 4% 

T = 1600 


a) < wo 

1.50 X 10“® 

315 X 10"‘ 


O 

3 

11 

3 

1.51 X 10"® 

369 X 10"‘ 



ence of the wavelength dependent photodiode response, 
because we monitored the laser power during wavelength 
sweeps with the same type of diode. The wavelength 
dependent discrepancy between the interference ampli¬ 
tudes estimated by our fits (see Figure indicates that 
we are limited by diffraction related effects. Models that 
include wavelength dependent mode propagation, as well 
as corrections due to cavity misalignments [21] can be 
applied in future studies to improve the understanding 
and accuracy of the reflection response. 

With the current models, at low frequencies and for 
the low-finesse cavity we can calculate 


Uc = \J (2u^o)^ + (^a)^ + (u 2„)2 -f (0.003Mfl)2 -I- {Uy_^Y. 

(33) 


Table IV shows the resulting combined standard un¬ 
certainties. We also include the corresponding values 
achieved for the high-finesse readout. With the cur¬ 
rent low-finesse device we reach uncertainties of 0.25%, 
which is similar to levels achieved at NMI calibration fa¬ 
cilities. The characterization of the optical response is 
currently limiting; applying more complex cavity models 
m can improve this in future studies. The uncertainties 
achieved with the high-finesse system are lower by almost 
an order of magnitude, however, they are still limited by 
the optical response characterization. One should note 
that the high-finesse cavity in Ref. m was also mode 
matched, using a curved mirror to form a hemispheri¬ 
cal resonator, potentially reducing the influence of beam 
propagation effects neglected in our Airy function analy¬ 
sis. 
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B. Bias stability 

For a given model of the dominating noise sources one 
can compute a corresponding Allan deviation behavior 
|29H33j ■ By comparing the two most dominant noise 
types, white noise and 1/f noise, one can compute the 
Allan deviation minimum and derive from that the ex¬ 
pected bias stability. 

A white noise of Hw = a™ • I/'s/Hz induces an Allan 
deviation slope of 


a„(r) = ^. (34) 

A 1/f noise of Oi/y = ai/f/f ■ 1/-\/Hz induces an Allan 
deviation slope of 

cri//(r) = 2TT°^yfr. (35) 

For the minimum value, which is often denoted as the 
bias stability, we can now calculate the optimal integra¬ 
tion time. 


"n/S CLuj 

Tbs = • 

ZTT ai/f 


(36) 


This leads to an estimate of the bias stability of 


^bs 


47r 




(37) 


Assuming a white acceleration noise floor of Ou, = 
8p.g/-\/Hz and a 1/f noise floor of a^, = 0.9mg///-\/Hz 
the integration minimum can be found at Tbs ~ 2.5 ms 
and the corresponding bias stability is on the order of 
0.2 mg, which is consistent with our data presented in 
Figure Using the presented formulas we can now ex¬ 
trapolate the 1/f laser frequency stability that is re¬ 
quired to achieve a certain bias stability with our model. 
This is calculated as 


A// 


dvr 


(38) 


With our cavity parameters and the earlier assumed 
white noise floor we require a 1/f laser frequency noise 
of less than /i/^///-\/Hz « 1700 Hz/f/-\/Hz to achieve a 
bias stability of 1 pg at integration times of 1.3 seconds. 


VI. SUMMARY AND CONCLUSION 

We have presented an accelerometer device that is 
a promising candidate for future applications requiring 


high bandwidth and in-situ self-calibration. We have con¬ 
ducted a detailed analysis of the device physics and de¬ 
rived the effects that influence noise, accuracy and stabil¬ 
ity. We achieve a bias stability of better than 0.3 x 10“^ g 
and an estimated scale factor uncertainty of 0.15%. For 
an input signal of 2 g we can estimate a total acceleration 
measurement error on the order of 3.3 mg, dominated by 
systematic effects in the characterization of the cavity 
response. 

The bias stability of our device is currently limited 
by laser frequency noise and non-stationary thermal and 
acoustic noise coupling into the measurement through 
parasitic beams in the fiber set-up. At frequencies above 
1 kFlz our readout noise is dominated by shot noise. The 
accuracy of our device is limited by the characterization 
of the reflectivity response of the Fabry-Perot cavity. 

The limits of the achievable accuracies for scale fac¬ 
tor and bias stability have been discussed, both for low- 
finesse and for high-finesse readout. While the ultimate 
limits of self-calibration to SI standards are rather well 
understood based on the applied simple behavior models, 
targeted studies will be conducted in the future to mea¬ 
sure the experimental variations of the scale factor over 
long time scales and its dependency on temperature. Ap¬ 
plying a cavity model that takes diffraction effects and 
misalignments into account may further improve the self¬ 
calibration of the cavity response. Further experiments 
are also required to test the limits of the simple behav¬ 
ior models by, for example, performing detailed measure¬ 
ments of the harmonic oscillator response function. 


Disclaimer 


Certain commercial equipment, instruments, or mate¬ 
rials are identified in this paper in order to specify the 
test and measurement procedure adequately. Such iden¬ 
tification is not intended to imply recommendation or 
endorsement by the National Institute of Standards and 
Technology, nor is it intended to imply that the materials 
or equipment identified are necessarily the best available 
for the purpose. 
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